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Abstract 

We study the homogenization of a steady diffusion equation in a highly 
heterogeneous medium made of two subregions separated by a periodic 
barrier through which the flow is proportional to the jump of the tem- 
perature by a layer conductance of the same order of magnitude of the 
materials in consideration. The macroscopic governing equations and the 
effective conductivity of the homogenized model are obtained by means of 
the two scale convergence technique. We show that under some hypothesis 
the homogenized systems contain convective terms of order one. 



1 Introduction 

Homogenization in niulticomponcnt e— periodic media with interfacial barriers 
has been extensively studied these last years. There are many mathematical 
works devoted to the subject and we refer the reader for instance to Auriault 
et al [5], [6], Benveniste [8], Canon et al. [3, Hummel [11], Monsurro [Hj and 
Pernin [T^ . . . 

In [5] , [Sj the layer conductance or sometimes called the resistivity is consid- 
ered as a positive function of order of magnitude and five distinct macroscopic 
models are derived by the formal asymptotic expansion method [7], [18j . These 
homogenized models are related to five values of 7 which are —2, —1, 0, 1 and 2. 
Monsurro [14] considered the same problem and the derivation of the homoge- 
nized models is obtained with the help of the oscillating test functions method 
of Tartar [T^]. The case 7 < — 1 has been studied by Hummel [llj but for me- 
dia with disconnected components arranged in a tesselation configuration. He 
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used the two-scale convergence method by Nguetseng [T2] , further developed by 
Allaire [3j and adapted to periodic surfaces, see for e.g. Allaire et al [4] and 
Neuss-Radu [TS] . The aim of this paper is to consider in the simplest configura- 
tion the case of resistivity having zero average value on the periodic interface. 
As in Ainouz [2] we show that the homogenized problem contains convective 
terms. 

The paper is organized as follows: in Section 2 we give the geometry of 
the medium in which the stationary diffusion micro-model is set. In Section 
3 the weak formulation of the problem is studied in a suitable function space 
and we give the main a priori estimates. Finally in Section 4, we recall the 
two-scale convergence and its main results, which we shall use it to derive the 
homogenized problem with the help of two-scale convergence method. 

2 Setting of the Problem and statement of the 
main result 

We begin this section by describing the geometry of the medium underlying 
the exact micro-model for the steady diffusion equation in a two-component 
periodic medium. 

Let Y = (0, 1)" be the unit cell of periodicity and assume that Y is divided as 
F = Yi UE where Yi,i = 1,2 are two open subsets of Y and E = 9Yi n9l2 is 
a sufficiently smooth interface that separates them. The sets Yi and I2 are made 
of two different materials but having conductance of same order of magnitude 
and E is a thin layer of very low conductance which constitutes a flow exchange 
barrier. 

Let Ai{y) (i = 1,2) denotes the conductivity tensor of the material Yi. We 
assume that Ai(y) is a n x n Y-periodic matrix- valued function and continuous 
on M" such that 

m,\X\^ <A{y)X- \ (2.1) 

and 



A,{y)X-fj<A'h\XM (2.2) 

for a.e. y ^ Yi and for all X,rj ^ M", where m^, Mj are positive constants. 

Let a denotes the barrier resistivity of E. For simplicity we suppose that a 
is a continuous function on M", Y— periodic with a zero average- value on E, i.e. 

a[y)da{v) = 0. (2.3) 

Here da{y) denotes the surface measure on E. We decompose a into its positive 
and negative parts as follows: 

a = a~^-a~, (2.4) 
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where 

= sup {a, 0) , a~ = sup {—a, 0) . 
We assume that 

a+ > ao (2.5) 

where ao is a positive real number. 

Let tti (y) {i = 1,2) be continuous functions on M" and F— periodic such 
that 

Oi (y) >Vi>0 for a.e. y e Yi. (2.6) 

Let f2 be a bounded domain in K" with sufficiently smooth boundary F. 

Let £ > be a positive number taking its value in a sequence of real numbers 
converging to zero. We consider the following open subsets of fl 

ni = {x€fl; Xi(J) =1} i = l,2 (2.7) 

where Xi iv) is the F— periodic characteristic function of Y^. 
The (n — l)-dimensional surface defined by 

= {a; e O; 3"fc e Z", ^ + ~k ^ ^} (2.8) 

is then by construction the interfacial barrier that separates the materials ilf 
and f2|, i.e. dfll n 9^1 . 

Let rf = Hf n r (i = 1,2) and assume that |rf| ^ 0. Let ff € ^^(Of) 
{i = 1, 2) be uniformly bounded functions, that is, 

llfflkoj <C (2.9) 

where (and throughout this paper) C is a positive constant independent of e. 
Let us consider the following transmission problem: 

-diY{M^)Vut{x))+ai{^)uUx)=f!{x) inQt, ^ = 1,2 (2.10) 
Ai{^)Vul{x) ■ = A^{^)Slu\(x) ■ on S^ (2.11) 



Ai{^)Vu\{x) ■ = -a(^) {u\ - ul) on E^ (2.12) 

< = Oonrf. i = l,2 (2.13) 

where is the unit outward normal to ilf obtained in an obvious way by 
extending y, the unit outward normal to Yi, by F— periodicity to the whole 
space M". 
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The first equation (|2.10p is the classical steady diffusion problem. The second 
equation (|2.11[) expresses the continuity of the flow across the interfacial contact 

while the boundary condition (|2.12p says that the flow is proportional to the 
jump of the temperature. For a physical justification of the condition (|2.12p we 
refer the reader to Carslaw and Jaeger [lOI, Kholodovskii [12j. Finally the last 
condition (|2.13p is the well-known Dirichlet boundary condition. 

In order to state the main theorem of this paper we begin by setting some 
notations. Let us first consider the following microscopic problems: 



-dWy{My){Vyj,{y)) = in Y,, (2.14) 

iMy)i^yl^iy)) ■ Ky) = -aiy) on S, (2.15) 

y I — > li{y) F-periodic, « = 1, 2, (2-16) 
and for k — 1, 2, . . . , n the problem 

-div,(A,(y)(e'= + ^yiKy)) = in Y,, (2.17) 

(A,(2;)(e'= + V,ef (y)) • ^{y) = on E, (2.18) 

y ^ (y) r-periodic, z = 1, 2, . (2.19) 



where div^ (resp.Vy ) is the divergence (resp. gradient) operator with respect 
to the variable y and {e'^)i<k<n is the canonical basis of R". 

The well-posedness of these elliptic problems is well-known, see for instance 
Bensoussan et al. [7]. Let us mention that the assumption (|2.3p is a compati- 
bility condition for the unique solvability of the problems (I2.14p - (|2.16p . 

Let us define the effective matrices A^^^ = {a^'')i<k,j<n and the convection 
vectors Bi ^ {bf)i<k<n (« = 1,2) 

= ^ Me" + V,ef ) • {e' + ^v^dy, (2.20) 

bf^i-iy-' [ A,Vj,dy+ [ a^^daiy). (2.21) 
Let us denote 

d - ^ «(7i - l2)da{y) (2.22) 

and 

c,=d + J aidy, gi{x) ^ J fi{x,y)dy, i = 1,2 (2.23) 

where fi is the two-scale limit of the sequence (xi(f )/i^)e>o- 
The main result of the paper is the following theorem 
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Theorem 2.1 Let {ul,ul) be the solution of the problem S2.10\) - i2.13\) . Then 
up to a subsequence, (xi(f X2(f )w2)e>o two-scale converges to (ui,W2) G 
(i?Q(f2))^ solution of the Homogenized Problem: 

-div{Al^-^\/ui) + Bi ■ Vmi - B2 ■ Vu2 + ciUi - du2 = gi in ft, (2.24) 

-div(^2''^-''VM2) + B2 ■ Vm2 - Bi ■ Vui + C2U2 - dui = 32 in ^, (2.25) 
ui = U2 = on r. (2.26) 

3 Solvability of the problem and a priori esti- 
mates 

Let us introduce the Hilbert space V = {^{,^1) x {^21^2) where 

H^iniTl)^{v,eH'ini);v,=OonTl}. i ^ 1,2 (3.1) 
is equipped with the norm: 

11(^1, «2)||y. := Ikilllnj + \\v2\\lni + l^i " ^^2|laE- (3-2) 

The equivalent weak formulation of the micro- model (|2.10p - (|2.13p is : 

For each e > 0, find (uf , ul) G such that , , 

a%iulu'2),ivuV2)) = L%{vi,V2)) for all ivi,V2) € V ^'^''^^ 

where for all {wi,W2), ivi,V2) € 

amwi,W2),{vi,V2)) = iA{-)Vw' -Vv^ + a{-)w''v^)dx + 

Jn £ £ 

+ / a{-){wi-W2)ivi-V2)da'{x), (3.4) 



Lmvi,V2)) = ^J^ rix)v'{x)dxj . (3.5) 
Here da'^{x) denotes the surface measure on and A, w^, u*^, f^ are defined by 
A{y) - xiiy)Aiiy) + X2(y)^2(y), y e M", 



= Xli-)wi{x) + X2{-)w2ix), X e n, 

X X 
V" = Xl{-)viix)+X2{-)v2ix), xen, 

e e 
r - xi{^)f^{x)+X2{^)fl{x),xen. 



To study the solvability of the problem (13. 3p we shall use the following inequality: 
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Lemma 3.1 There exists a constant Ci > independent of e such that for 
every Vi G (Jlf ) and for all Si > we have 

<CK;^lb.||^,oj+'5»e||Vz;,||2,,j). *-l,2 (3.6) 

Proof. See for instance Ainouz [2] (see also Monsurro 14 ). Note that the 
constants Ci depends solely on the geometry of F^. ■ 
Throughout this paper we assume that 



lin^^-i^. (3.7) 

•1=1,2 Oi 

where Ci are the constants defined in Lemma 13.11 Now we give the following 
existence result : 

Proposition 3.2 Let the assumption Jg. 7| ) be fulfilled. Then for any fixed e > 
0, there exists a unique couple (itf,?/!) € solution of i3.3\) . Moreover we 
have the a priori estimate 

WiuluDWv. <C. (3.8) 

Proof. First we show that a^(-, •) is coercive on . Let {vi, V2) G V^, then 

a'{ivi,V2),ivi,V2)) > Ej=i,2(™d|V«»llo,nf +'?dl^^'llo,nf) ^3 
+ ("0 - ||q:"||oo,eO /s. (■"1(2:) - V2{x)f d(j^{x). 

But in view of Lemma 13.11 we see that 



Set 



Therefore (|3.9p becomes 

a'{ivi,V2),{vuV2)) > 51 ('^■'ll^^'llo.Jif ++^''ll^*llo,f2f) 

+ao\\vi - V2\\l^SE . (3.10) 

Now we choose, for example, Si = — - — ^^^'^ i = 1,2. Then by 

r]imi + (Ci\\a ||oo,eO 

(|3.7p we have 

^ _ m^jrj.m, - (C,||q~||oo,sO^) ^ ^ 
77imi + (Gi||a ||oo,sO 



/ (viix) - V2{x)f da%x) < J2 ai^Mll^. +S.e\\Vv,\\ln.). 



i=l,2 



C 
Oi£ 
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and 



^2 



Vimi-{Ci\\a ||oo,s^, , . 

2Cimi\\a ||oo,E= 

Hence by p. lip and p.l2p . the inequality (|3.10p becomes 

a"((wi,f2),(vi,i'2)) > co||(wi,W2)||y. (3.13) 

where 

Co = min(/3i,/32,7i,72,Q;o) > 

which is independent of e. 

It is easy to observe that a'^{{'Wi, W2), ('^i, ^'2)) is bilinear continuous on x 
V' and that L^((wi,W2)) is linear and continuous on V^. Consequently by the 
Lax-Milgram Lemma the problem p.3p has a unique solution (lif,?/!) in . 
Furthermore we have 

\\iul,ul)\\h < -L^iulu^)) < 1 ^ ll/f ||o,oj IKIIcOf 

which implies by (j2.9p that the sequence (wf,U2)£>o is uniformly bounded in 
. The Proposition is then proved. ■ 
In view of the a priori estimate p.Sp . one is interested in investigating the 

limit of the sequence ((uf,it|))e>o as e ^ in a sense that will be specified 

later. This is the purpose of the next Section. 



4 two-scale convergence process 

In this Section we shall use the two-scale convergence method to determine the 
Homogenization of the Problem (|2.10p - (|2.13p . For more details on two-scale 
convergence, we refer the reader to Nguesteng [TB], Allaire [3] and a recent 
paper by Lukkassen et al. [13] ■ We also mention the work by Allaire et al. 
[1] (see also Ainouz [T], Neuss-radu [TS]) where the two-scale convergence is 
applied to periodic surfaces. For convenience we recall the definition and the 
main compactness results of this method. 

Definition 4.1 

1. Let (w^)e>o be a sequence in L'^{fl) and vq e L^(n x Y). We say that 
iv^)s>o two-scale converges to vq if for every ip S I?(r2;C^r (Y)) we have 

lim / v"^ {x)(p{x, —)dx = / vo{x,y)(p{x,y)dxdy. (4.1) 
^^°Jn £ JnxY 
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2. Let {v'^)^yQ be a sequence in L^(I]'^) and vq £ L^(ri X E). We say that 
{v^)e>o two-scale converges to vq if for every Lp G T){Q]C^j. (Y)) we have 

lime / v'^{x)(p{x,—)da'^{x)= / vo{x,y)ip{x,y)dxda{y). (4.2) 

We now give the main compactness results of the two scale convergence : 
Theorem 4.2 

1. Let {v'^)s>o be a uniformly bounded sequence in L^(f7). Then one can 
extract a subsequence which two-scale converges to a function vq G L'^{flx 
Y) in the sense of Definition \4 ■ j^ l 

2. Let {\^v^)e>o be a uniformly bounded sequence in L'^{Yf). Then one 
can extract a subsequence which two-scale converges to a function wq € 
L^{n X S). 

3. If {v^)e>o is uniformly bounded in H^{fl) (resp. _ffg(f2)j then one can 
extract a subsequence still denoted {v^)e>Q md there exist v € H^(fl) 
(resp. H^in)) and V G L^{n; H^^^(Y)/R) such that 



i) (w^)e>o two-scale converges to v in the sense of Definition [^7 

ii) (Vw^)£>o two-scale converges to + V^V in the sense of Definition 
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) Moreover {\/£v^)i,yo is uniformly bounded in L^(S'^) and {v'^)e>o two- 
scale converges to v in the sense of Definition\4.^ we have wq = vq. 



Next we shall apply these results to determine the two-scale limit of (u?)e>o 
and the source terms (/f )j>o, i = 1, 2. 

Lemma 4.3 For each i = 1,2 there exist fi <E L'^{Q x Yi), Ui{x) G Ho{^) 
and Ui G L^(0, i/pg^(li)/R) such that up to a subsequence we have for all 
if, G P(l^;Cp- (Yi)) and G V{n;C^^^ (F,))"; 



lim / fi{x)ipi{x,^)dx ^ J^^J^ fi{x,y)lp^{x,y)dydx, (4.3) 



e->0 



lini / u^{x)(pi{x,^)dx = J J Ui{x)(pi{x,y)dydx, (4.4) 



lim / \/u^{x)^i{x,-)dx ^ / / [\/ui{x) -\-\/yUi{x,y)]iJi{x,y)dydx (4.5) 
and for all ip G V{Ti;C^^ (Y)) 
lim / eul{x)ip{x, —)da^ (x) — / / Ui{x)ip{x,y)da{y)dx, i — 1,2. (4.6) 
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Proof. The existence of the two-scale limits is an immediate consequence 
of the a priori estimates (|2.9[) . (|3.8p and the definition of two-scale conver- 
gence. Indeed, since for i = 1,2 the sequences Xi(f)/f(a;), Xi{j)uf{x)s>o^ 
{Xi{j)^u1{x))s^o and {^/euKx))s>o are uniformly bounded in L^(Sl), L'^{Q), 
(L2(0))" and ^^(se) respectively. Then by Theorem US one can extract a 
subsequence still denoted e and there exist fi{x,y), u^{x,y) € i^(ri x Y) and 

G X Y)Y such that Xi(f )/f (a^), Xii'Mi^) and x.(f )Vuf(x) two- 

scale converge respectively to fi{x,y), Ui{x,y) and £,i{x,y). We point out that 
fi{x,y), Ui{x,y) and S,i{x,y) are equal zero outside F^. Arguing as in Allaire 
[21 Thm 2.9] we easilty arrive at u^{x,y) = where Ui{x) e Hq{^) 

and C»(a;,y) = VMi(x) Vj,i7,(a;, y) where U, € L^{n, H^^^(Y,)/R). Finally, by 
Allaire et al. [1 Prop. 2.6] (see also Ainouz [1] or Neuss-Radu [15, Thm. 2.2])) 
we obtain the last limit. ■ 

The following result is an extension of an auxiliary result given in 

Lemma 4.4 Up to the subsequence given in Lemma \4.3\ we have 
lim / uf{x)ipt{x,-)da^{x) = / / Ui{x,y)ifi{x,y)da{y)dx 



for any ipi £ ^{Q-jCper {Yi}) such that Jj^ipi{x,y)d(T{y) =0 i = 1,2. 

Proof. For a fixed x in fl, and i = 1,2 let (j)i{x,y) be the solution of the 
following boundary value problem 

-divy(l),{x,y) =0 ii-iYi (4.7) 

(j>i{x, y) ■ v{y) = (-1)'"^ ipt{x, y) on E, (4.8) 

y I — > (j>i{x,y) is F-periodic. (4-9) 

Such a function exists since Jj^ipi{x,y)da{y) — 0. Furthermore, the solution 
(j)i{x,-) belongs to Hp^j.{Yi)/R. Taking into account the boundary conditions 
and (|4.9p we see that 

f X f X 

ul{x)Lpi{x,y)da^{x) ^ I ^ul{x)4>i{x,-)dx+ I uf{x)diva;{<t'iix,-))dx. 

(4.10) 

By (|4J| and letting e -> in (|4T0| together with (|4?5l) we have 

lim / ul{x)ifi{x,y)da%x) = [\/Ui{x) +\/yUi{x,y)](l)i{x,y)dydx 



Ui{x)divx<j)i{x,y)dydx. (4.11) 
'n JYi 
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Since Ui G H}^(p.), integration by parts with respect to x gives 
\Iui{x)(t)^{x,y)dx + I Ui{x)divx(t)i{x,y)dx ^ 



Therefore (|4.1ip becomes 

hni / u1{x)(p,{x,y)da''{x) ^ I I VyUi{x,y)(t>i{x,y)dydx (4.12) 
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Again, integrating by parts with respect to y in the right hand side of (|4.12p 
and using ((48)) . ((49)) yield 



(-1) / U^{x,y){(j)^{x,y) ■ i^{y))da{y)dx 
Ui{x,y)ipi{x,y)da{y)dx. (4-13) 



VyU^{x,y)(j)i{x,y)dydx = - Ui{x,y)diYy(f),{x,y)dydx + 

Yi JYi 



Finally, by combining ()4.12p and (|4.13p we arrive at the desired result. This 
proves the Lemma. ■ 
Now we are in a position to give the two-scale homogenized problem. 

Proposition 4.5 The two-scale limit {ui,Ui) G i?o(fi) x L^{Q, Hp^j.{Yi)/M.) is 
the solution of the two-scale homogenized system 

{ ~ d:iY y {Ai(y)('sj u,{x) + y yU dx, y))) = mVl-KY,, i=l,2 
-div^(/^^ Aj(VMi + VyUi)dy) + /y, ai{y)ui{y)dy+ 
+ i-iy-' a{Ui - U2))da{y) ^g^mU, z = 1, 2 
(^i(Vui + VyC/i)) • V = (^2(Vu2 + VyC/2)) -v onVL^H, 
(Ai(Vui + VyC/i)) • V = -a{ui - U2) = on n X T,, 

Ui — on T , y 1-^ Ui{x, y) is Y -periodic. i — 1,2 
Proof. Let (p^{x) e V{n), $i e V{n,C^^{Yi)) i ^ 1,2. We take Vi{x) 

X 

. + e^i{x, — ) in the weak formulation (|3.3p . We have 

e 

a'iiulu'2),{vi,V2)) = /[ + /| + /| + /| + eJf + eJ| 
Lmvi,V2)) = Kl + KI + sLl+eLl 

where 

(a,(-)V< (a;)(V(/3,(x) + Vy<I>j(x, -) + a,{-)u^^{x)ipi{x)\ dx, i 
I e e e J 
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(4.14) 



1,2 



Jt.^ e e e 

Jf= / A,{-)Vul{x)V^^r{x,-)dx, i = l,2 

11= I ot{-){u\{x) - ul{x)){ipi{x) - ip2{x))d(T'' {x); 



and 



Kt= / f!{x)^^{x)dx■, Ll= / f!{x)<S>,{x,-)dx, i = l,2. 

Letting e ^ and using (|4.3p - (l4.6p . we have 

limZf = J^^Y^ {A,{y){^u,{x) + \/yU^{x,y)){\/(p^{x) + \/y<^>^{x, y))dxdy} + 
InxY, ai{y)ui{x)(pi{x)dxdy, i = 1,2 

(4.15) 



lim/| = / a{y){ui{x)-U2{x)){^i{x,y)-^2{x,y))dxd<7{y){A.W) 
limi^f = / f,ix,y)^^{xdxdy- 1 = 1,2 (4.17) 



By virtue of Lemma 14.41 

lim/|=/ a(tj)(C/i(.T,?;)-C/2(a;,?/))((pi(a;)-(^2(a;))da;dCT(?/); (4.18) 



On the other hand we see that 



I Jf I + |L^| < C, z = l,2 

where C is a positive constant independent of e. Hence 

hmejf = hmeif = 0, i = 1,2 (4.19) 

Now passing to the hmit in (|3.3p and using (|4.15p - (|4.19p . we obtain the two-scale 
system 

V / {A,{Vui + \7yU.i){\7ipi + Vy^i) + a^u,ipi\ dxdy+ 



1=1,2 



{a{Ui - U2){(pi - (P2) + a{ui - U2)($i - $2)} dxda{y) 



/ fiLpidxdy 



(4.20) 
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Now by density the formulation (|4.20p remains true for all (fi, W2, Vl, V2) in 
{H^{n))^ X L^{n,H^^^{Y,)/m) x L^{n,H^^^{Y2)/R). An integration by parts 
yields the two-scale homogenized system (|4.14l) . ■ 

Proof of Theorem [HH First we take $i,(^2,*2 = in (|42Q1). Then 
integration by parts yields 

-diVa; JY^Ai{y){Vui{x) +VyUi{x,y)dy + Jy^ai{y)ui{x)dy 

+ J^a{y)iUi{x,y) - U2{x,y)da{y) = /^^ fi{x,y)dy in n, (4-21) 
ui = on r. 

Similarly we take tpi,^i, $2 equal to zero and this gives 



— diVa; 



Jy, ^i(y)(V'"i(2;) +^yUi{x,y)dy + J ai{y)ui{x)dy 



+ I^(^iy)iUi{x,y) - U2{x,y)da{y) = /^^ fi{x,y)dy in fJ, (4 22) 

lii — on r. 



Now choosing ipi,ip2,^2 —0 (resp. ipi,ip2,^i = 0) in (|4.20p gives after integra- 
tion by parts 



~dWy{Ai{y){Vui{x) + VyUi{x,y)) = in f7 x F^, 

(Ai(y)(Vwi(a;) + ^yUi(x,y))) ■ v^{y) 

= {A2{y){\/U2{x) + VyU2{x, y))) ■ vi{y) onflxT., 

{Ai{y){Wui{x) + VyUi{x,y)) ■ vi{y) + a{y){ui{x) - U2{x)) = on f7 x S, 

y\ — > Ui{x,y),U2{x,y) F-periodic. 

(4.23) 



Next we shall decouple the problem (|4.23p , that is eliminating the unknowns 
Ui,U2 from the system (I4.14p . The linearity of the problem and the fact that 
Ui do not depend on the fast variable y enable us to put 

C/,(x,y) = ^X.^(y)^(a^)+7.(2/)(ui(a;)-M2(a;)) + M.(a:), i = l,2 (4.24) 



k=l 



Then inserting (jiTM)) into (g^Bland KT]\ yields the equations 
Thus we have proved the Theorem 12.11 

Remark 4.6 It is easy to see that if Ai are symmetric for all i then Bi ~ 0. 
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